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Abstract 

Purpose: The traditional brachytherapy catheter reconstruction with biplane images is based on digitizing radio- 
opaque markers with a pointing device on a film or on a screen. An algorithm to automate digitization of radio-opaque 
marker coordinates on biplane images is presented. 

Material and methods: To obtain the marker coordinates in a proper sequence, instead of usual pair of recon- 
struction images, series of images were taken with insertion of radio-opaque markers consecutively into the catheters. 
The images were pre-processed to suppress the shield of anatomic structures. The determination of the marker coordi- 
nates is based on the detection of characteristic high gradient variation in pre-processed image profiles. The method was 
tested with six endometrial insertions performed with Simon-Norman catheters using our version of Heyman packing. 

Results: 28 catheters of six treatment fractions were digitized, typically 10 markers per catheter. To obtain the mark- 
er coordinates, adjustment of two threshold levels on the pre-processed images were needed. The coordinates of the 
radio-opaque markers on the biplane projection images were obtained without positive or negative artefact. 

Conclusions: The dummy source coordinates on the biplane images were digitized in a proper sequence: from 
the catheters' tip towards the end of the catheters. After the three-dimensional reconstruction of the catheters from 
the digitized coordinates, the geometry file was imported by the brachytherapy planning system for dose calculation. 
The method has the advantage to eliminate manual digitization of the dummy sources. 
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Purpose 

Multiple flexible brachytherapy catheters are applied 
by inserting the breast with endometrial treatments, using 
modified Heyman packing method with inflatable rectal 
applicator, and with surface insertions. Reconstruction 
of large number of catheters with biplane projection im- 
ages is a time consuming process. More methods were 
published to speed up treatment planning with decreasing 
number of digitized points, and with elimination of the 
need of one-two-one coordinate correspondence require- 
ment [1]. The use of automated applicator reconstruction 
has not become a common practice with the commercial 
brachytherapy planning systems. Only two methods were 
published in this field: one of them by Li et al. for recon- 
struction of Fletcher-Suit applicator [2], and the other by 
Rosenwald for localization of curved wires [3]. 

The modified Heyman packing according to Herbol- 
sheimer et al., usually is performed with 6-11 pieces of 
Simon-Norman catheters. The heads of the catheters are 
uniformly distributed in the uterine cavity [4]. At our in- 
stitution, endometrial insertions with smaller amount of 



Simon-Norman catheters (Nucletron, an Elekta company, 
Elekta AB, Stockholm, Sweden) are performed. Anoth- 
er significant factor is the heads of each Simon-Norman 
catheters are positioned to the fundus uteri [5] as illus- 
trated in Figure 1. Depending on the dimension of the 
uterus, the number of catheters is usually between 5 or 6. 
The dose distributions are optimized with dwell times to 
fit the uterine contour [6]. The geometrical arrangement 
of the actual insertion approximates the ideal one illus- 
trated in Figure 1. It differs more or less from that, there- 
fore we prepare individual treatment plan for each treat- 
ment fraction. Even with reduced number of catheters, 
the number of points to digitize is over 50. An algorithm 
to automate the digitization of radio-opaque markers in- 
serted into the Simon-Norman catheters is presented. 

Material and methods 

Obtaining biplane reconstruction images 

At our institution, the reconstruction images are ob- 
tained with image capture from the video signal of a mo- 
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bile C-arm image intensifier (Giraldoni, Italy). With the 
small power C-arm, better image quality was experienc- 
ed using PO images between 150 degrees and 142.5 de- 
grees, than with usual lateral ones [7]. The fluoroscopy 
image size was 512 x 512 pixels with 0.32 x 0.32 mm 
resolution, the pixel values were represented with 8 bits. 
The images were inverted, therefore the markers are 
appeared with high pixel values. The Simon-Norman 
catheters inserted into the uterus are close together and 
we experience overlap on the reconstruction images. 
The traditional biplane imaging was not suitable for the 
digitizing algorithm. To obtain the marker coordinates 
in a proper sequence (first catheter: first marker... n lh 
marker, second catheter: first marker... n th marker, etc.), 
we have taken series of PA and PO images with markers 
inserted consecutively into one of the catheters. One of 
the images is illustrated in Figure 2A. 

Digitization of the marker coordinates 
on the reconstruction images 

Before digitizing the dummy sources, we suppress 
the pixel values of anatomical structures as follows. We 
obtain two copies of each reconstruction image: one of 
them is translated laterally by +1 pixel, while the other 
by -1 pixel, and the two images are subtracted. On the 
reconstruction image, the neighbouring pixel values is 
similar, therefore the anatomical structures on the differ- 
ential image are suppressed (see Figure 2B). In contrary, 
the profiles taken across the markers show high gradi- 
ent variation as illustrated in Figure 2C. The maxima and 
minima of the profile at the marker vary with shielding 
of bony structures, and with decreasing sensitivity of 
the image intensifier at the image circumference. On the 
differential images, a low and a high threshold levels are 



Fig. 1. Our version of Heyman packing: ideal arrangement 
of five Simon-Norman catheters with 6 mm diameter heads, 
each catheter is inserted to the fundus uteri 



adjusted until each radio-opaque marker appears on the 
image without artefacts (see Figures 2B and D) coordi- 
nates of the markers, obtained with an in-house written 
scanner program, which reads the differential image 
rows consecutively and searches for minima and max- 
ima ('M' and 'm' in Figure 2C) of the profiles between 
the sections determined by the two threshold levels. 
The location of a radio-opaque marker on an image row 
is the mean value of position of minima and maxima. 
Depending on the length of the radio-opaque marker, 
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Fig. 2. Posterior-oblique reconstruction image of radio-opaque markers inserted into a Simon-Norman catheter (A), the dif- 
ference image with a low and a high (L and H) threshold levels (B), the pixel values smaller than the low threshold level are 
represented with black, while pixel values larger than the high one are showed with white. Profile obtained across a marker 
along the section shown on the differential image (C) and histogram of the pixel values, and the threshold levels (D) 
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the position is detected on 2-15 consecutive profiles. 
The coordinates of a radio-opaque marker are computed 
by averaging the numbers of the consecutive rows, where 
the marker is detected and by averaging the correspond- 
ing lateral position values. The image manipulations 
are performed with the ImageMagick software package 
(http://www.imagemagick.org), and the differential im- 
ages with the actual threshold levels are visualized with 
the ImageJ program http://www.rsbweb.nih.gov/ ij). 

Determination of the optimal translation 

To determine the optimal value of translation of the 
reconstruction images, we performed the image process- 
ing with shifts by 1 to 5 pixels. The standard deviation 
of the pixel values in a region out of the markers (N), 
the difference (S) of the maxima (M) and minima (m) of 
the profiles taken across a selected marker, and the ratio 
S/N on the differential images were computed and illus- 
trated in Figure 3. 

The influence of the catheter inclination 

The Simon-Norman applicators follow the shape of 
the uterine cavity, therefore part of the catheters are not 
perpendicular to the image profiles. We simulated the 
catheter inclination with rotating the reconstruction im- 
age of a strait catheter through 15, 30, and 45 degrees. 
The number of marker detections in consecutive image 
profiles was counted and summarized in Table 1. 

The influence of the image quality 

The quality of the reconstruction images taken with 
the corpulent patient is often poor. We investigated 
the influence of the image noise on the digitizing process. 
We obtained noisy images from a good quality recon- 
struction image with adding noise with the Image} pro- 




gram illustrated in Figure 4. The S/N values on the cor- 
responding differential images (A...D) were: 34, 28, 12.8, 
and 10, respectively. 

Testing the method with Simon-Norman 
applicators 

The marker digitizing method was tested in insertions 
with 3-6 pieces of Simon-Norman catheters of 6 mm dia- 
meter head, as illustrated in Figure 1. Since the original 
markers supplied with the Simon-Norman catheter set 
were not suitable for the method, we replaced them with 
radio-opaque markers 1 cm apart. We digitized 28 cathe- 
ters of 6 insertions, and 10 markers per catheters. We com- 
puted the quantities S, S/N and the difference of the high 
and low threshold levels (H-L). These quantities were aver- 
aged for the PA and the PO images separately summarized 
in Table 2. 

Results 

The quantities S and N of the differential images in- 
creased with the value of the translation, while the ratio 
S/N decreased. The optimal value of translation was two 
pixels. Inclination of the catheter through 45 degrees de- 
creased the number of marker detections in consecutive 
profiles by less than 3; the corresponding uncertainty was 
±1 pixel. The lateral uncertainty of the marker positions 
of the non-inclined catheter was also ±1 pixel. The S/N 
values of the images in Figure 4 were 28.7, 24.7, 12.8, and 
10, respectively. The markers were recognized without 
artefact on the images A, B and C in Figure 4, while there 
were a missed marker and more artefacts on the very noisy 
image D. In insertions performed with the Simon-Norman 
catheters, the PA reconstruction images showed good or 
medium quality (S/N > 60 and S/N > 30, respectively), 
while the quality of the PO images was poorer. All mark- 
ers were recognised without positive or negative artefact, 
the short markers in 2-4 consecutive profiles, while the 
long ones in 10-14 profiles. The difference between the 
high and low threshold levels H-L were 30 times the value 
of N for good quality images, which decreased to 20 and 
10 for medium, and poor ones. 

Table 1. The number of marker detections in conse- 
cutive image profiles with rotated catheter 



Fig. 3. The signal amplitude (S) and the standard devia- 
tion of fluctuation (N), and S/N of the differential image 
versus the translation of the reconstruction image 
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Fig. 4. Reconstruction images with added noise 



Discussion 

Weitman et al. described the use of the modified Hey- 
man packing with MR images. They inserted 11 pieces 
of Simon-Norman catheters that were uniformly distri- 
buted in the uterus. Since the catheters identification 
was not possible with sectional images, the catheter re- 
construction was performed with biplane X-ray images. 
The reconstructed coordinates were combined with the 
contours obtained with the MR [9]. At our institution, the 
endometrial treatments are performed with 5 or 6 piec- 



es of catheters positioned to the fundus uteri, the most 
probable location of the endometrial cancer. In our pre- 
vious study, we showed that with this method higher 
dose can be delivered to the fundus uteri than with the 
Y shaped applicator [5]. With dwell time optimization, 
acceptable dose coverage can be achieved for the medium 
size uterus [6]. In lack of access to sectional imager, we 
apply a C-arm as a low cost alternative of the integrated 
brachytherapy unit. We capture the fluoroscopy images 
and convert to DICOM format, and store to a pictorial 
archive. Since the brachytherapy planning systems (BPS) 
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do not allows to import the dummy sources coordinates, 
therefore we reconstructed the three-dimensional coordi- 
nates of the catheters. However, analytical solution can 
be applied in this case and we used our multi-parametric 
fit method [8]. After reconstruction of each catheter, we 
generated a geometry file of 3D coordinates, which was 
imported by the Plato BPS vl3.7 (Nucletron, an Elekta 
company, Elekta AB, Stockholm, Sweden) for dose cal- 
culation. 

Conclusions 

The coordinates for the three-dimensional reconstruc- 
tion of the dummy sources should be obtained with 
a proper sequence on the reconstruction images. Since 
the Simon Norman catheters are close together and often 
overlapped, this requirement was ensured with taking 
a series of biplane images with markers inserted con- 
secutively into the catheters. The images were pre-pro- 
cessed and scanned with profiles close perpendicular 
to the catheters. The dummy source coordinates were 
obtained from the catheters' tip towards the end. Small 
intervention of the user to adjust the threshold levels on 
the pre-processed images was required. However, for the 
method a series of images, which number is two times the 
number of catheters is needed, has the advantage of elim- 
ination of manual digitization of large number of marker 
points and ensures accurate reconstruction of Simon Nor- 
man catheters. 

Acknowledgements 

The Heyman packing insertions for the study were per- 
formed by Kofi Agyemang-Prempeh M.D. 

Disclosure 

Author reports no conflict of interest. 
References 

1. Li S, Chen GTY, Pelizzari CA et al. A new source localization 
algorithm with no requirement of one-to-one Source corre- 
spondence between biplane radiographs. Med Phys 1996; 23: 
921-927. 

2. Li S, Pellizzari CA, Reft C et al. Computer-aided reconstruc- 
tion of Fletcher-Suit Source positions. Med Phys 1994; 21: 
1123-1130. 

3. Rosenwald JC. Automatic localization of curved wires used 
in brachy therapy. Comput Prog Biomed 1975; 4: 103-112. 

4. Herbolsheimer M, Sauer O. HDR brachytherapy for endome- 
trial cancer. International Nucletron Radiotherapy Journal 1996; 
Special Report No.8: 34-45. 

5. Palvolgyi J, Agyemang-Prempeh K. Modified Heyman pack- 
ing performed with small number of Simon-Norman cathe- 
ters. / Contemp Brachytherapy 2010; 2: 33-36. 

6. Palvolgyi J. Dose Distribution of Modified Heyman Packing. 
Physica Medica 2006; 22: 127-130. 

7. Palvolgyi J. To what extent can digital images obtained with 
a non-isocentric C-arm be used for brachytherapy treatment 
planning in gynaecology. Radiother Oncol 2006; 67: 107-112. 

8. Palvolgyi J. Multiparametric fit method in reconstruction of 
Fletcher-Suit-Delclos applicator. Med Phys 2006; 33: 69-75. 



9. Weitmann HD, Potter R, Waldhaus C et al. Treatment of en- 
dometrial carcinoma with 3D image-based high dose rate 
brachytherapy using modified Heyman packing: clinical 
experience and dose volume histogram analysis. Int } Radiat 
Oncol Biol Phys 2005; 62: 468-478. 



Journal of Contemporary Brachytherapy (2014/volume 6/number 1) 



